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Suppression of 3D-ordering by defects in the S' = 1/2 frustrated chain magnet LiCuVCb 
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We report on a heat capacity study of high quality single crystal samples of LiCuVCK - a frustrated 
spin S = 1/2 chain system - in magnetic field amounting to 3/4 of the saturation field. At low 
fields up to about 7 T, a linear temperature dependence of the specific heat, C p oc T, resulting from 
ID magnetic correlations in the spin chains is followed upon cooling by a sharp lambda anomaly of 
the transition into a 3D helical phase with C p oc T 3 behavior at low temperature. The transition 
from a spin liquid into a spin-modulated (SM) phase at higher fields occurs via a hump-like anomaly 
which, as the temperature decreases further turns into a C p oc T 2 law distinctive for a quasi-2D 
system. We suggest an explanation for how nonmagnetic defects in the Cu 2+ chains can suppress 
3D long-range ordering in the SM phase and leave it undisturbed in a helical phase. 

PACS numbers: 75.50.Ee, 76.60.-k, 75.10.Jm, 75.10.Pq 


Quantum-spin chains with frustrated exchange inter¬ 
actions were among the most interesting issues for both 
experimental and theoretical research in condensed mat¬ 
ter physics of the past decade (dj-iS]. The enhanced effect 
of quantum fluctuations imposed upon a fine balance of 
exchange interactions leads to a variety of novel ground 
states and phase transformations in these systems HE- 
LiCuVCK is an example of a quasi-ID magnet whose un¬ 
conventional magnetic phases result from a competition 
of ferromagnetic and antiferromagnetic exchange inter¬ 
actions between nearest-neighbor (Ji) and next-nearest 
neighbor (J 2 ) in-chain spins. As a result of this partic¬ 
ular combination of exchange interactions a helical in¬ 
commensurate spin structure is stabilized in this system 
below T)v « 2.3 K jg]. A strong reduction in the ordered 
spin component of Cu 2+ ions in this state {p)/hb — 0.3 
Si 0] provides evidence that the system partially retains 
properties of ID chains. 

Moderate applied magnetic fields of 7-8 T induce a 
transformation of the spin helix into a collincar spin- 
modulated (SM) structure in which all the spins are par¬ 
allel to the field, with their ordered components oscillat¬ 
ing along the chain with an incommensurate period.This 
transition may be related to the field evolution of short 
range chiral (transverse) and spin-density wave (longitu¬ 
dinal) correlations for the ID J 1 -J 2 model [13- 

In the field range just below the saturation field 
(HoH sa t — 45 T) 11| the theory predicts the presence 


of a long-range nematic ordering |7j, LL2|]. The magnetic 


properties of LiCuVCK in this field range were studied 
by pulse magnetization and nuclear magnetic resonance 
(NMR) techniques |ll|, ll3] • These experiments specified 
the supposed held range where the spin nematic phase 
can exist. The experiments discussed in our report were 


performed at fields well below this range. This work is 
solely focused on the magnetic properties of LiCuVCK in 
the helical and SM phase at lowest temperatures, and the 
spin-liquid phase for elevated temperatures. 

Single crystals of LiCuVCK grown by the methods de¬ 
scribed in [14( belong to the same series as the ones stud¬ 
ied in our previous high Held NMR experiments E3- 
The main part of the specific heat measurements were 
carried out in the Tallahassee National High Magnetic 
Field laboratory with the use of a home-made relaxation 
type calorimeter operating at temperatures 1.8-30 K in 
magnetic Held up to 35 T. Supplementary measurements 
were performed on a Quantum Design Physical Property 
Measurement System (PPMS) with a 3 He insert and a 
9 T cryomagnet. The samples of 0.67 and 5.8 mg by 
mass were studied in the first and the second experiment 
respectively. The magnetic field was applied along the 
c axis of the crystal. 

Zero-field results obtained in both experiments and 
presented in Fig. [T| are in perfect agreement with each 
other. The phonon contribution to the heat capacity (ex¬ 
pressed in units J/(K mol)) has been approximated by 
the function C p = 1.4 • 10~ 4 T 3 (solid line in the upper 
panel). It was chosen so that the magnetic entropy (with 
the phonon part subtracted) attains the value of 0.71? In 2 
at a temperature 30 K in agreement with numerical re¬ 
sults for the ID J 1 -J 2 model 0- This rough estimation 
gives a correction to the total heat capacity of about 3%, 
8% and 20% at T = 5, 10 and 15 K, respectively, shown 
by the error bars in the lower panel of Fig. [T] Within 
these bounds our data conform with previous measure¬ 
ments 16, 17]. This panel shows a comparison between 
the magnetic part of the specific heat C m (T) obtained 
at zero field and at a field of 9 T in the PPMS measure- 
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FIG. 1. (Color online) Upper panel: temperature dependence 
of the specific heat at zero magnetic held measured using 
a PPMS (circles) and using a home made high-held setup 
(squares); the solid line is an approximation of the phonon 
contribution to the heat capacity. Lower panel: magnetic 
specihc heat C m (T) (phonon contribution subtracted) of the 
sample measured in a PPMS at = 0 (O) and 9 T (•), 
the error bars are explained in the text, solid line is a linear 
ht to both curves in the temperature range from 5 to 10 K; 
the inset shows the specihc heat divided by temperature as a 
function of T obtained at several values of magnetic held (see 
legend). 


ments. No significant difference between the two curves 
is observed in the high-temperature range. All the data 
in the lower panel of Fig. [T] exhibit temperature depen¬ 
dences C m oc T in the range 5 — 10 K providing evidence 
for the presence of ID short-range correlations in the spin 
chains. The linear in temperature section of a zero-field 
curve is followed by a sharp lambda anomaly at the tran¬ 
sition into a 3D helical phase at Tjy = (2.27 ± 0.03) K in 
agreement with a number of previous studies [§, 9j. The 
inset to the lower panel of Fig. [T] presenting C m /T(T) 
curves at various fields traces the evolution of this transi¬ 
tion shifting to low temperature as the field increases. For 
the held hqH = 9 T which is greater than the held of the 
transition into a SM phase a lambda anomaly is replaced 
by a hump-like feature followed by a C m oc T 2 depen¬ 
dence towards lower temperatures. This curve slightly 
deviates from the square law at lowest temperatures. 

This low-temperature part of the data requires more 
careful investigation and is expanded in the upper panel 
of Fig. [2] The phase transition at zero held is followed 
by a distinct C m oc T 3 behavior of the specihc heat with 


FIG. 2. (Color online) Upper panel: Expansion of the 
low temperature part of the inset to Fig. 1, solid lines are 
quadratic and linear extrapolations to the zero field and 9 T 
curves respectively, dashed-dotted line is a calculated contri¬ 
bution from 63 Cu and 65 Cu nuclei at the field 10 T. Lower 
panel: integrated C m /T(T) curves at various field values nor¬ 
malized to 7?ln2; inset shows the boundary of the helical 
state obtained from both the PPMS and the high field mea¬ 
surements, the line is a guide-to-the-eye. 

a small linear contribution revealing itself at the lowest 
temperature (corresponding ht is presented in the Fig. [2j. 
Shifting of the Neel temperature under magnetic held 
(TV values are annotated in the upper panel and shown 
by open circles in the inset to lower panel) is accompa¬ 
nied by the increase of this low-T contribution. In the 
absence of the transition peak at [ 1 qH = 9 T the high- 
T linear part of the curve is converted into a quasi-2D 
C m oc T 2 law extending for half a decade by temperature 
(see linear ht to C m /T curve for 9 T). The low-T rise of 
this curve might indicate the proximity to a 3D tran¬ 
sition. Extrapolation of the phase boundary HqH C 2 (Tn) 
(see the inset to the lower panel of Fig. 0 to zero temper¬ 
ature gives the value /io-ffc2(0) -9T. Alternatively, this 
rise could be attributed to the specihc heat of the Cu nu¬ 
clei (the corresponding curve estimated for the external 
held 10 T and zero hyperhne held is shown in Fig.0. In 
any case, the total entropy evolved in the process of low- 
T disordering does not depend on the external held in 
this held range: the corresponding functions obtained by 
integration of C m /T(T) curves converge at temperatures 
above the ordering transition at zero held. 

A further increase in the applied magnetic held leads 
to a considerable change in the specihc heat behavior (see 
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FIG. 3. (Color online) Magnetic part of the specific heat di¬ 
vided by temperature as a function of T measured at various 
magnetic fields using a high field setup; shadowed area cor¬ 
responds to the formation of static 2D magnetic correlations 
as observed in an NMR experiment HI , vertical lines mark 
temperatures of scans in magnetic field shown in the inset. 
Solid line in the inset is a linear fit to data at T = 5.2 K in 
the intermediate field range. 

Fig. [3]). A broad hump in C m /T vs T curve preceding 
the formation of static 2D correlations shifts to some¬ 
what higher temperatures, while the quasi-ID linear part 
(constant for C m /T) breaks down at fields above 20 T. 
It is replaced by another broad feature merging with the 
first one at a maximum experimental field of 33 T. Two 
scans in an applied magnetic field performed at constant 
temperature below (1.55 K) and above (5.2 K) the or¬ 
dering transition T)v = 2.27 K at (H = 0) are presented 
in the inset to Fig. [3] The low-T curve has a sharp peak 
associated with the transition from a helical phase at 
^qH C 2 = 6.0 T followed by a gradual decrease in mag¬ 
nitude towards higher fields. The specific heat reaches 
its minimum around 20 T and then starts to increase as 
the field enters the range of critical fluctuations near sat¬ 
uration. In contrast, the high-T curve demonstrates a 
crossover from a nearly constant to a linear field depen¬ 
dence C/T(H) (the fit is shown by solid line) extending 
to about 25 T. The final downturn of the specific heat 
developing at fields above 30 T indicates the approach of 
the system to the saturated state. 

We shall now discuss the specific heat data obtained for 
high-quality single crystal samples of LiCuVCH at field 
and temperatures of spin-liquid and ordered helical and 
SM phases. The specific heat in the spin-liquid phase for 
temperatures 5 < T < 10 K is a nearly linear function 
of T at low fields (H < H sat / 4) and becomes nonlinear 
at higher fields. The observed evolution is in qualitative 
agreement with numerical results obtained for a Heisen¬ 
berg S = 1/2 chain model jTI]. The temperature depen¬ 
dence of the specific heat in the ordered helical state can 
be roughly approximated as a sum of linear and C m oc T 3 
terms while for the SM state the cubic term is replaced 
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FIG. 4. Schematic diagram of the influence of defects on 
the spiral state (a) and the spin-modulated state at fioH ~ 
20 T (b), defects are shown by black dots. The spiral state 
is nearly insensitive to the chain break while the phase of a 
spin-modulated chain is randomly shifted at the point of the 
defect. 

by Cm oc T 2 contribution. Different exponents in power 
laws point to 3D and quasi 2D properties of these phases 
respectively confirming that the antiferromagnetic order¬ 
ing between neighboring a6-planes in the SM structure is 
destroyed. This effect was previously reported in neutron 
scattering and NMR measurements |20j,[2]j. The replace¬ 
ment of a sharp lambda anomaly at the transition to the 
helical state by a broad feature at the crossover to the 
SM state gives another hint of this change. 

In the following we propose a qualitative model to ex¬ 
plain this different dimensional behavior by considering 
the effect of a nonmagnetic impurity. A composition 
study of LiCuVCH reveals a Li deficiency of a few per¬ 
cent even in crystals grown using the low temperature 
technique which provides crystals of the highest perfec¬ 
tion. It was argued that the Li deficiency results in holes 
with spins S' = 1/2 localized on oxygen which in turn 
form a Zhang-Rice singlet with the neighboring copper 
spin [3, El. Such a singlet should be equivalent to a 
non-magnetic impurity replacing a Cu spin and will pro¬ 
voke an unusual magnetic state in its vicinity with the 
two nearest ferromagnetic exchange bonds removed, but 
with a conserved antiferromagnetic bond between the two 
parts of the interrupted chain. The presence of such de¬ 
fects in LiCuVCH single crystals of the same series as in 
this work was suggested in order to explain the bulk mag¬ 
netization process, as well as the vanadium NMR spectra 
in the vicinity of the saturation field [ijj}]. A schematic 
diagram of helical and SM spin chains broken by a non¬ 
magnetic defect is shown in Fig. [4] In contrast to a spin 
chain with only a nearest-neighbor exchange interaction, 
the frustrated J\-Ji chain is not fully broken by a non¬ 
magnetic impurity. The magnetic moments on both sides 
of the defect are coupled by indirect antiferromagnetic ex¬ 
change interaction leading to a local distortion and phase 
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shift of an incommensurate exchange structure. 

Obviously, the presence of this type of defect does 
not significantly disturb the helical magnetic structure 
in LiCuV 04 since the incommensurate wave vector of 
the helix k; c = 0.468 • 27 x/b (the unit cell of LiCuV 04 
contains two magnetic ions of Cu) corresponds to the 
angle between next-nearest neighbor spins close to ir (see 
Fig. QJi) and does not depend on the applied field. A 
residual phase shift of a helical structure near the defect 
may be eliminated due to an intra-chain exchange inter¬ 
action so that the initial 3D long range ordering is not 
suppressed. In contrast, the influence of defects on the 
SM ordering is expected to be more pronounced for the 
following reasons. The wave vector of a SM structure in 
a ID J 1 -J 2 model is field dependent and expressed as fol¬ 
lows: k; c = (1 — M{H)/M sa t) -7T /b 0,13-Q. This relation 
has been traced experimentally in the field range below 
15 T [23,0] and therefore, the value ki c is significantly 
less than 7 x/b at all fields above H c 2 . The elongation of 
the spin modulation period along the chain in field leads 
to a random shift of its phase at each defect (see Fig. HJd, 
the scheme of the moment distribution over the chain is 
presented for fioH ~ 20 T). The boundary conditions for 
spins at both ends of a finite chain should depend on its 
length leading to an uncertainty in the phase shift of an 
incommensurate structure due to a random distribution 
of the vacancies. As in case of a helical structure, the 
effect of these inter-chain exchange interactions is to re¬ 
store these broken correlations and establish a 3D order¬ 
ing. Taking into account the hierarchy of exchange pa¬ 
rameters (the interaction between c planes is the smallest 
one) correlations between spins in neighboring ab planes 
are hindered and the system remains in an effectively 2D 
ordered state. 

In conclusion, a specific heat experiment was per¬ 
formed on high-quality single crystals of LiCuV 04 . The 
results show that while a helical spin state observed in 
the low-field range H < H c 2 clearly demonstrates 3D 
magnetic ordering, the field induced magnetic phase at 
H Z > H c 2 has properties peculiar to a quasi-2D system. 
We propose a simple model to describe how nonmag¬ 
netic defects in the Cu 2+ chains can destroy 3D long- 
range ordering in a SM phase, but leave it undisturbed 
in the helical phase. In addition, the proposed model 
allows one to simultaneously explain the results of both 
NMR 0, 0, 241 and neutron scattering measurements 


0 resolving the existing discrepancy in their interpreta¬ 
tions. According to NMR data, the magnetic phases at 
both sides of H c 2 are dipolar so that the average val¬ 
ues of magnetic moments of Cu ions are non-zero on 
a time scale of NMR measurements (~ 1 sec). This 
proves the presence of static magnetic correlations and 
precludes an identification of the spin-nematic state, as¬ 
sumed in 0 for higher magnetic fields H > H c 2 . Nev¬ 
ertheless, the short-range character of these correlation 
especially in the direction perpendicular to ab planes de¬ 


tected in neutron diffraction patterns can be understood 
using the above model describing the suppression of a 
3D long range order by non-magnetic defects. Further 
specific heat experiments extending to the vicinity of the 
saturation field (about 45 T) is especially interesting in 
prospect of possible observation of the transition from 
a SM into a theoretically predicted quadrupole nematic 
state. 
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